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We discuss the structure of a central stellar cluster whose dynamics is influenced by gravitation of 
a supermassive black hole and by the dissipative interaction of orbiting stars with an accretion disc. 
We also take the effect of disc self-gravity into account. We show that the cluster properties are 
determined predominantly by the radial profile of the disc surface density. To this aim we develop 
a simple steady-state model of the central cluster and we estimate the rate at which stars migrate 
to the centre. 

This model is relevant for central regions of an active galactic nucleus (AGN) containing a rather 
dense accretion medium. In passing we also briefly mention a possibility that a fossil (inactive) disc 
could exist in the centre of our own Galaxy. Such a hypothetical disc could perturb the motion 
of stars and set them on highly elliptic trajectories with small pericentre distances. The required 
mass of the disc is less than one percent of the central black hole mass, i.e. below the upper limit 
permitted by present accuracy of the orbital parameters of S-stars in the Galactic Centre. 



I. INTRODUCTION 

We examined the long-term orbital evolution of 
stars forming a dense stellar cluster surrounding a cen- 
tral black hole with an embedded accretion disc. This 
configuration is relevant for central regions of active 
galactic nuclei 0, 0| and it may be applied also to 
the center of our own Galaxy, assuming that rapid ac- 
cretion took place and a gaseous disc was formed at 
some stage of its history 8]. Indeed, the presence of a 
dissipative gaseous environment can provide a mech- 
anism driving stars towards the black hole, while the 
gravitational influence of the disc may pump eccen- 
tricity of the orbits to large values at some moments 
of the evolution This would help us to under- 

stand the puzzling nature of young stars in the close 
vicinity of the black hole in the centre of our Galaxy. 
An attractive feature of this calculation is that it pro- 
vides a well-defined model allowing to estimate the 
expected time-scales of the orbital migration as well 
as the distribution of eccentricities. It may turn out 
to be more likely that a different (non-standard) kind 
of a gaseous disc or a dusty torus plays the role, but 
this would not change the essence of the model. 



II. RESULTS 

We idealise a galactic core as a system consisting of 
a central black hole, an accretion disc and a dense stel- 
lar cluster. The three components interact with each 
other. Naturally, various approximations need to be 
imposed in order to keep our model sufficiently simple 
and tractable. The aim is to examine the structure of 
a stellar system in the region of black hole gravita- 



tional dominance i?h, including the effects due to a 
gaseous disc. Two regions of the cluster can be dis- 
tinguished according to the characteristic time-scale 
of processes dominating the stellar motion. The outer 
cluster is assumed to reach a gravitationally relaxed 
form [H, acting as a reservoir of fresh satellites that 
are being continuously injected inwards. The inner 
cluster is defined as a region where star-disc collisions 
take over. 



Individual orbits 

Stars loose their orbital energy and momentum by 
means of successive dissipative passages through an 
accretion disc. If the orbit is inclined with respect 
to the disc plane, each encounter with the disc slab 
slightly modifies the orbital parameters. The over- 
all trend is to circularize orbits and to decline them 
into the plane of the disc 0. Characteristic time 
of ahgning stellar orbits with the disc is td(a) ~ 
to Ms {Th./Tiq) ^ (a/i?g)' yr , where a is semi-major 
axis, to and q are constants (order of unity for the 
standard disc model) and Ms = M,/{10^Mq) is the 
fiducial value of the central mass in AGN 4] . 

Different modes of radial migration apply to stars 
in the disc plane, depending whether a star succeeds 
to open a gap in the disc medium, or if it remains 
embedded entirely and proceeds via density waves ex- 
citation. We switch between relevant modes in the 
numerical integration. 

A gap is cleared in the disc if the embedded star is 
sufficiently massive and its Roche radius exceeds the 
characteristic vertical thickness of the disc slab at a 
corresponding radius, « (M*/M)^/'^r > h{R). In 
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FIG. 1: An example of the episode during which the steUar orbit becomes highly elliptical. Left: occasional jumps toward 
high eccentricity occur due to the Kozai-type mechanism in the field of the central black hole and a gravitating disc. 
Prevailing circularisation is expected because of continuing energy losses and the resulting orbital decay via star-disc 
collisions, however, fraction of highly eccentric orbits always persists in the cluster, as indicated by the colour scale. 
These trajectories bring stars close to the black hole where they can be preferentially captured or destroyed. Right: the 
ratio of mean latitudinal frequency to the mean radial epicyclic frequency uir (for the same orbit as on the left). It 
turns out that the jump of orbital parameters occurs at the moment when the ratio u>z/uJr = 1. The indicated time-scale 
corresponds to the case of M, = 4 x IO^Mq and a hypothetical disc with the mass ~ O.OIM, . 



this case the motion of the star is coupled with the 
disc inflow, 

flgap = /(a,AiE,M,;&,gj) (1) 

with 

/(a, fiE, M.; b, q,) = ^bMl' y^^^a^- (^^^ 

with /ZE being the accretion rate in units of the Ed- 
dington accretion rate and a viscosity parameter. The 
numerical factor b and power-law indices (71...4 are de- 
termined by details of the particular model adopted 
for the disc 

On the other hand, if the star is too tiny to create 
the gap, than density-waves are the dominating mi- 
gration process. The resulting orbital decay can be 
written in the form 

fldw = 77^ /(a,ME, (2) 

Mq 

Differences in the rate of stellar migration are thus 
introduced in the model already within this very sim- 
plified picture where the process of satellite sinking is 
driven by the gas medium. The dependence on the or- 
bital parameters, the stellar masses and sizes of stars 



causes the gradual segregation of different stellar types 
in the nuclear cluster. 

The gravitational field of the disc provides a per- 
turbation capable of exciting large fiuctuations of ec- 
centricity. Assuming that these fluctuations occur on 
the time-scale substantially longer than the orbital pe- 
riod, one can apply an averaging techique and evolve 
equation for mean orbital parameters. See e.g. Kozai 
I3 and Lidov & Ziglin for a general introduction 
to the formalism; see Vokrouhlicky & Karas 14] and 
Subr & Karas JJJ for its modification to the present 
situation. We computed the gravitational field of the 
disc and took its effect into account for the long-term 
orbital decay of stellar trajectories. The effect of ec- 
centricity oscillations is shown in Figure ^ Here, we 
have set the current (time zero) values of orbital eccen- 
tricity e{t) and semimajor axis a{t) identical as those 
reported for S2 star in Sagittarius A*. For the central 
mass we adopted A'/, = 4x lO^Af0. A false-colour plot 
shows what fraction of time the orbiting star spends 
with a certain value of the orbital eccentricity. There- 
fore, this graph illustrates the expected fluctuations of 
eccentricity and explains the existence of highly ellip- 
tic trajectories, assuming that an accretion disc or a 
torus was present at the centre of our Galaxy. The re- 
quired mass of the disc Md was of the order of fraction 
of percent of the central black hole mass, consistent 
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FIG. 2: Left panel: distribution of semi-major axes throughout the stellar cluster, modified by interaction with the gas- 
pressure dominated standard accretion disc. The broken power-law profile is established by different types of interaction 
governing different regions of phase space. The three stellar sub-samples are present in our model and they can be easily 
distinguished in the plot: the outer reservoir (magenta), the inner part of the cluster which is dragged by the disc (blue), 
and a subsample of stars fully embedded in the disc (green). Right panel: Graph of mean inclination (i) in the cluster. 
The reservoir is spherically symmetric ((i) = 7r/2), the dragged cluster is somewhat fiattened (0 < (i) < 7r/2), and the 
embedded population is located in the disc plane ((i) — 0). 



5^ -1 



CO 
O 




0.005 0.010 
Line-of-sight velocity Iv^l/c 



0.015 




0.005 0.010 
Line-of-sight velocity \v^\/c 



0.015 



FIG. 3: Velocity profile along the line of sight of the inner cluster, integrated across a column of cross-sectional radius 
Rb = 10''7?g (left panel) and Rb — lO^iig (right panel). Blue/red lines represent different view angles of the observer: 
^ = 0° and 60° , respectively. Growing anisotropy of the modified cluster produces the dependence of measured line profile 
on 5, i.e. 1=1 (vz',^)- The green line is for the reference Bahcall-Wolf distribution, which exhibits spherical symmetry. 
See ref. [Tl| for further details. 



with the present-day upper limit. Characteristic time- 
scale tc of the oscillations depends on M^/M, ratio, 
and is indeed much longer that the dynamical time at 
corresponding radius (of the order of thousand orbital 
revolutions in our case, Md/M, < 0.01). 



The cluster 



Now we construct a steady-state cluster assuming 
that it is supplied with fresh stars from the reser- 
voir at a rate inversely proportional to the relaxation 
time. The effects introduced in the previous section 
are taken into account. Our computational scheme 
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FIG. 4: The spatial density n«(r) and the corresponding projected density of the cluster are shown using logarithmically 
spaced levels of shading. Columns (a) and (b) represent the meridional section at two different scales, namely, 4Rd and 
0.27?d across (radii are expressed in terms of the disc outer radius, -Rd). Next columns are the edge-on (c) and the 
face-on (d) projections of the cluster. Across columns, the upper row shows the referential cluster P) (n* cx r"'^''*). In 
subsequent rows, the system has been already modified via the interaction with two types of discs, case (iii) and case 
(ii)-l-(iv), as discussed in ref. pl|. 



allows us to further distinguish between two subsam- 
ples of the inner cluster: the dragged inner cluster 
consists of stars on orbits crossing the disc periodi- 
cally; the embedded inner cluster is formed by stars 
entirely aligned with the disc. See Subr et al. [ij for 
details of our approach. For the definiteness of exam- 
ples we assumed that the central mass (M, = lO^M©) 
is surrounded by the gas-pressure dominated Shakura- 
Sunyaev disc (s — —3/4) with M = O.lMEdd and 
viscosity parameter a — 0.1 (so these values are ade- 
quate for an AGN). The outer stellar cluster can be 
characterized by the number density uq = lO^pc"'^ 
and velocity dispersion CTc = 200km/s. 



Structure of the cluster modified by the 
interaction with an accretion disc 



Figure|21shows the density structure of the modified 
cluster. Majority of stars forming the embedded clus- 
ter sink to the centre in the regime of density waves, 
hence Vr oc r^/^ and n(a) oc a~^/^. In the dragged 
cluster, the orbital decay leads to the governing in- 
dex given by s — 1/2 = —5/4, and the corresponding 
number density n(a) oc a^^^. The asymptotic profile 
oc a^^^ of the outer cluster is determined by the initial 
distribution. Isotropy of the outer cluster is violated 
in the inner regions where the mean inclination satu- 
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rates at « 0.2 (dragged cluster). The influence of the 
disc manifests itself in different characteristics of the 
inner cluster. For example, dependence of the drag 
on size and mass of stars causes gradual segregation 
of different stellar types present in the cluster. 

FigureOlshows the integrated properties of the clus- 
ter that can be compared with observation. We plot 
the shape of a synthetic spectral line I{vz), i.e. inten- 
sity in the line as a function of line-of-sight velocity Vz 
near the projected center of the cluster. Local maxi- 
mum of the line occurs around v ~ sin(^)z;K(i?d). For 
some values of model parameters, this secondary peak 
exceeds the central maximum and dominates the pre- 
dicted profile. High-velocity tails of the line profiles 
are also noticeably affected in comparison with an un- 
perturbed form of the outer cluster pj| . 



Oblateness of the cluster 

Figure 21 shows various two-dimensional sections of 
the cluster arranged in four columns. One can clearly 
observe the impact that star-disc collisions have on 
the cluster structure, namely, an increasing oblateness 
of the stellar population in the core and, in some cases, 
the tendency to form an annulus of stars. The reason 
for different structures is the continuous crashing of 
stars on the disc plane. Furthermore, in case of stars 
embedded in the disc, different modes of star-disc in- 
teraction occur and facilitate their radial transport to 
the central hole. 

It is worth noticing that the gravity of the oblate 
cluster adds with the field of the embedded accretion 
disc. This way the effect of orbital oscillations, dis- 
cussed in previous sections, may be further enhanced. 
Indeed, any non-spherical perturbation, for example 
general-relativity effects of frame-dragging near a ro- 
tating black hole, can potentially be responsible for 
similar kind of oscillations. However, the issue cannot 
be resolved within the framework of our simple model, 
because other effects intervene (for example the frag- 
mentation of the disc into selfgravitating clumps) and 
these may act in the opposite direction. 



III. DISCUSSION 

The orbits are aligned and circularized at typical 
radii of 10^ 10^i?g, spiralling further on nearly cir- 
cular orbits towards the centre. This provides limits 
on gravitational waves emerging from the cluster. The 
rate of the capture events can be estimated as 



The total accretion rate onto the black hole is a sum 
of Ms (which involves massive stars in the cluster), 
and the accretion rate M (gas in the disc). 

It is worth noticing that the orbital decay of stars 
near a black hole is indeed relevant for forthcom- 
ing gravitational wave experiments, because the gas- 
dynamical drag should be taken into account with suf- 
ficient accuracy. It was estimated 0,01 that this effect 
can be safely ignored at late stages, shortly before the 
star plunges into the hole, if accretion takes place in 
the mode of a very diluted fiow, but the situation is 
quite different in case of AGN hosting rather dense 
nuclear discs. Gas-dynamical effects most likely dom- 
inate over the gravitational radiation as far as their 
influence on the cluster structure is concerned. 

The average rate of energy losses which the orbiter 
experiences via gravitational radiation (per revolu- 
tion) can be written in terms of orbital parameters 
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7/2 



/i(e). 



(4) 



Linked with this £ is the change of semi-major axis 

-3 



h = -1.28 X 10^'^cM. 



M, 







hie) (5) 



and the loss of angular momentum 
32 Af 5/2M2 



{!I2 (1 _ g2y 



/3(e). 



(6) 



Functions /i , fi and /a are of the order of unity. The 
above-given formulae Q-® assume the star is on an 
elliptic trajectory in Schwarzschild geometry. It there- 
fore provides only an order-of-magnitude estimation, 
because collective effects operating in the cluster are 
neglected. 

Gravitational radiation orbital decay competes with 
the orbital decay caused by hydrodynamics of star- 
disc encounters. The relative importance of these two 
influences depends on the ratio between the hydro- 
dynamical drag (which is roughly proportional to the 
disc density) versus gravitational-wave losses (which 
increase with orbital eccentricity). Discussing this in- 
terplay is a tricky task once the gravitationally in- 
duced oscillations of eccentricity are taken into ac- 
count. See refs. 0,01 for a more detailed discussion of 
these effects. Substantially different results and time- 
scales can be expected for objects with rather dense 
discs (such as the case of standard discs in AGNs) and 
those which are inactive and represent a rather clean 
system (such as our own Galactic Center). 



Ms 



lO^^M 



5/4 



no 



106 pC- 



Mr, 



V 104i?g 



Mq yr- 



(3) 



IV. CONCLUSIONS 

Although the above-described model is intended 
mainly for active galactic nuclei with a relatively dense 
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accretion disc or a dusty torus, the problem of stars 
crashing on a gaseous disc may be relevant also for 
the centre of our Galaxy. The modified cluster struc- 
ture is relevant for estimating the rate of black-hole 
feeding and, vice-versa, for the issue of feedback that 
a super-massive black hole exhibits on the the host 
galaxy. 

The inherent limitations persist in our present dis- 
cussion, namely, we have not incorporated a fully self- 
consistent treatment of the disc gravity. Even though 
we computed the gravitational field across a suffi- 
ciently large domain of space, we did not account for 



the feedback which star-disc collisions exert on the 
disc structure. Introducing some kind of a clumpy 
model of the disc will be very interesting, as it may 
exert a more substantial effect on the cluster structure 
by elevating the impact of star-disc collisions. For fur- 
ther details we refer the reader to papers [ill IS i 
and to references cited therein. 
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